1 


High C/O Chemistry and Weak Thermal Inversion in the 
Extremely Irradiated Atmosphere of Exoplanet WASP- 
12b 

Nikku Madhusudhan 1 . Joseph Harrington 2 , Kevin B. Stevenson 2 . Sarah Nymeyer 2 , 
Christopher J. Campo 2 , Peter J. Wheatley’. Drake Deming 4 , Jasmina Bleeic 2 , Ryan A. 
Hardy 2 , Nate B. Lust ". David R. Anderson'. Andrew Collier-Cameron 6 . Christopher B. 
T. Britt 2 , William C, Bowman 2 , Leslie Hebty. Coel Hellier 4 , Pierre F. L. Maxted 4 , Don 
Pollacco 8 , Richard G. West 9 

Kavli Institute for Astrophysics and Space Research , and Department of Earth, 
Atmospheric , and Planetary Sciences, Massachusetts Institute of Technology, 
Cambridge , MA 02139, USA. Planetary Sciences Group, Department of Physics, 
University of Central Florida, Orlando. FL 32816-2385, USA. s Department of Physics, 
University of Warwick , Coventry, CV4 7AL.UK. '' NASA's Goddard Space Flight Center, 
Greenbelt, MD 20771-0001, USA. ' Astrophysics Group, Keele University, Staffordshire 
STS 5BG, UK. 6 School of Physics and Astronomy, University of St. Andrews, North 
Haugh, Fife KYI 6 9SS, UK Department of Physics and Astronomy, Vanderbilt 
University, Nashville, TN 37235, USA. ’ Astrophysics Research Centre, School of 
Mathematics & Physics, Queen's University, University Road, Belfast , BT7 INN, UK. 

9 Department of Physics and Astronomy, University of Leicester, Leicester, LEI 7RH, 
UK 

The earbon-to-oxygen ratio (C/O) in a planet provides critical information about 
its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 
causes a carbide-dominated interior as opposed to the silicate-dominated 
composition as found on Earth 1 ; the solar C/O is 0.54 2 . Theory shows that high 
C/O leads to a diversity of carbon-rich planets that can have very’ different 
interiors and atmospheres from those in the solar system 13 . Here we report the 
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detection of C/O a 1 hi a planetary atmosphere. The transiting hot Jupiter WASP- 
12b 4 has a dayside atmosphere depleted in water vapour and enhanced in methane 
by over two orders of magnitude compared to a solar-abundance chemical 
equilibrium model at the expected temperatures. The observed concentrations of 
the prominent molecules'* 6 CO, CHy, and H 2 0 are consistent with theoretical 
expectations for an atmosphere with the observed C/O = 1. The C/O ratios are not 
known for giant planets in the solar system, although they are expected to equal 
the solar value’ 8 . If high C/O ratios are common, then extrasolar planets are likely 
very different in interior composition, and formed very differently, from 
expectations based on solar composition 1,3,8 , potentially explaining the large 
diversity in observed radii. We also find that the extremely irradiated atmosphere 
(> 2500 K) of WASP-12b lacks a prominent thermal inversion, or a stratosphere, 
and has very efficient day-night energy circulation. The absence of a strong 
thermal inversion is in stark contrast to theoretical predictions for the most highly 
irradiated hot-Jupiter atmospheres 4 10,n . 

The transiting hot Jupiter WASP- 12b orbits a star slightly hotter than the Sun (6300 K) 
in a circular orbit at a distance of only 0.023 AU, making it one of the hottest exoplanets 
known . Thermal emission from the dayside atmosphere of WASP- 12b has been 
reported using the Spitzer Space Telescope 12 , at 3.6, 4.5, 5.8, and 8 pm° wavelengths, 
and from ground-based observations in the J (1.2 pm), H (1.6 pm), and K.s (2.1 pm) 
bands 14 (Figure 1). 

The observations provide constraints on the dayside atmospheric composition and 
thermal structure, based on the dominant opacity source in each bandpass. The J, H, and 
Ks channels ! 4 do not have significant molecular absorption features, and hence probe 
the deepest layers of the observable atmosphere, at pressure (P) of - 1 bar, where the 
temperature (T) is ~ 3000 K (Figure 1). The Spitzer observations 13 , on the other hand. 



are excellent probes of molecular composition. CH 4 has strong absorption features in 
the 3,6 pm and 8 pm channels, CO has strong absorption in the 4.5 pm channel, and 
H 2 O has its strongest feature in the 5.8 pm channel and weaker features in the 3.6 pm, 
4.5 pm, and 8 pm channels. The low brightness temperatures in the 3.6 pm (2700 K) 
and 4.5 pm (2500 K) channels, therefore, clearly suggest strong absorption due to CH 4 
and CO, respectively. The high brightness temperature in the 5.8 pm channel, on the 
other hand, indicates low absorption due to H 2 O. The strong CO absorption in the 4.5 
pm channel also indicates temperature decreasing with altitude, since a thermal 
inversion would cause emission features of CO in the same channel with a significantly 
higher flux than at 3.6 pm <,,lh . 

The broadband observations allow us to infer the chemical composition and 
temperature structure of the dayside atmosphere of WASP-1 2b, using a statistical 
retrieval technique*’. We combined a 1-D atmosphere model with a Markov-chain 
Monte Carlo sampler J that computes over 4x10° models to explore the parameter 
space. The phase space included thermal profiles with and without inversions, and 
equilibrium and non-equilibrium chemistry' over a wide range of atomic abundances. 
Our models include the dominant sources of infrared opacity in the temperature regime 
of WASP- 12b 5,18,19 : H 2 O, CO, CH 4 , CO 2 , lh~ H 2 collision induced absorption, and TiO 
and VO where the temperatures are high enough for them to exist in gas phase 9 ' 20 . The 
host star has a significantly enhanced metailicity (2 x solar) 4 , and evolutionary 
processes can further enhance the abundances s ; Jupiter has 3 x solar C/H (Ref. 7). Our 
models therefore explore wide abundance ranges: ~ 0.05 - 100 x solar for C/H and O/H, 
and 0.1 - 10 x solar for C/O. Figure 2 show's the mixing ratios of H 2 0, CO, CH 4 , and 
CO 2 , and the ratios of C/H, O/H and C/O, required by the models at different levels of 
fit. Figure 3 presents the temperature profiles. 
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We find a surprising lack of water and overabundance of methane (Figure 2), At 
2000 - 3000 K, assuming solar abundances yields CO and H 2 0 as the dominant species 
besides H? and He 19 ' 20 . Most of the carbon, and the same amount of oxygen, are present 
in CO, and some carbon exists as CFU. The remaining oxygen in a hydrogen-dominated 
atmosphere is mostly in H 2 0; small amounts are also present in species such as CCf. 
The CO/H 2 and H 2 O/H 2 mixing ratios should each be > 5 x 10' 4 , CH 4 /H 2 should be < 
10' , and CO 2 /H 2 should be - 10' , under equilibrium conditions at a nominal pressure 
of 0.1 bar. The requirement of H 2 O/H 2 s 8 x 10 ' 6 and CH 4 /H 2 a 4 x ! 0' 6 (both at 3a 
significance: Figure 2) is therefore inconsistent with equilibrium chemistry' using solar 
abundances. 

The observations place a strict constraint on the C/O ratio. We detect a C/O > 1 at 
3a significance (Figure 2). Our results rule out a solar C/O of 0.54 at 4.40. Our 
calculations of equilibrium chemistry' with a C/O ratio of 1 yield mixing ratios of H 2 O, 
CO and CH 4 that are consistent with the observed constraints. We find that, for C/O = 1 , 
FbO mixing ratios as low' as 10 ' ; and CH 4 mixing ratios as high as 10 " 5 can be attained 
in the 0.1 ~ 1 bar level for temperatures around 2000 K and higher. And, while the CO 
mixing ratio is predicted to be > 1 O' 4 , making it the dominant molecule after H 2 and Fie, 
CO 2 is predicted to be negligible (<30 w ). These theoretical predictions for a C/O = 3 
atmosphere, are consistent with the observed constraints on ffO, CH 4 , CO, and CO 2 
(Figure 2). 

The observations rule out a strong thermal inversion deeper than 0.01 bar (Figure 
3). Thermal inversions at lower pressures have opacities too low to induce features in 
the emission spectrum that are resolvable with current instruments. For comparison, all 
stratospheric inversions in solar system giant planets, and those consistent with hot 
Jupiter observations, exist at pressures between 0.01 - 1 bar 646 ’" 1 . The major 
contributions to all the observations come from the lower layers of the atmosphere, P > 



0.01 bar, where we rule out a thermal inversion (Figure 1 of SI). T he observations also 
suggest very efficient dav-night energy redistribution (Figure 2). The low brightness 
temperatures at 3.6 and 4.5 um Imply that only part of the incident stellar energy is 
reradiated from the davside, while up to 45% is absorbed and redistributed to the 
nightside. The possibility of a deep thermal inversion and inefficient redistribution was 
suggested recently 14 , based on observations in the J, H, and Ks channels, but the Spitzer 
observations rule out both conditions. 

The lack of a prominent thermal inversion contrasts existing theories, which 
designate WASP- 12b to the class of very hot Jupiters which are expected to host 
inversions 9,22 . At T > 2000K, molecules such as TiO and VO, which are strong 
absorbers in the UV/visible, are expected to be gas phase and potentially cause thermal 
inversions 9 . WASP- 12b therefore presents a major challenge to existing atmospheric 
classification schemes for exoplanets based on thermal inversions 9 ' 22 . However, it 
remains to be seen if the high continuum UV flux expected for WASP- 12b might be 
efficient in photo-dissociating inversion-causing compounds in the atmosphere thus 
explaining the lack of a strong inversion 10 . Alternatively, the present result might 
support theories which suggest that TiO and VO are less likely to be present aloft in the 
atmosphere to cause thermal inversions 20 . 

If high C/O ratios are common, then the formation processes and compositions of 
extrasolar planets are likely very different from expectations based on solar system 
planets. The host star has enhanced metallicity but die C/O ratio is unknown 4 23 24 . In the 
core accretion model, favoured for the formation of Jupiter, icy planetesimals 
containing heavy elements coalesce to form the core, followed by gas accretion 8,23 . The 
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abundances of elemental oxygen and carbon are enhanced equally'’ . maintaining a C/O 
like the star's. The host star WASP-12 is thus expected to have a C/O ~ 1. If the stellar 
C/O is independently determined to be < 1, then the C/O enhancement in WASP-12b ! s 
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atmosphere would suggest either an unusual origin for the planetesimals, a local over- 
density of carbonaceous grains'’' 2 * 1 , or a different formation mechanism entirely. 
Although carbon-rich giant planets like WASP- 12b have not been studied, theoretical 
studies predict myriad compositions for carbon-dominated solid planets 1 "’. Terrestrial- 
sized carbon planets, for instance, could be dominated by graphite or diamond interiors, 
as opposed to the silicate composition of Earth 1 " 1 . In the context of hot Jupiters, a 
dominance of carbon among the heavy elements in the interiors would likely change the 
mass and radius estimates from what have been calculated, based on solar abundances. 

The observed molecular abundances in the dayside atmosphere of WASP- 1 2b 
motivate a new' regime in atmospheric chemistry. It remains to be seen if 
photochemistry in WASP- 12b can significantly alter the composition in the lower layers 
of the atmosphere, P ~ 0.1 - 1 bar, which contribute most to the observed spectrum 
(Figure 1 of SI). Explaining the observed composition as a result of photochemistry 
with solar abundances would still be challenging. CH 4 is more readily photo-dissociated 
compared to H 2 O 11 " 2 '’, and hence a depletion of CH4 over that predicted with solar 
abundances might be expected, as opposed to the observed enhancement of CH 4 . Apart 
from the spectroscopically dominant molecules considered in this work, other species 
such as OH, C 2 H 2 , and FeH (Refs. 27, 28), which are not detectable by current 
observations, could potentially be detected with high-resolution spectra in the future. 
Most models of exoplanetary atmospheres have typically assumed solar abundances 
and/or solar C/O, thereby exploring a very' limited region of parameter space 926 ' 29 . The 
ease of WASP- 12b is a strong motivation for models to depart from solar abundances 
and abundance ratios. 
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Figure 1: Observations and model spectra for dayside thermal emission of 
WASP-12b. The black filled circles with error bars show the data: four Spitzer 
observations 13 (3.6 urn, 4.5 urn, 5.8 um, and 8 urn; see Table 1), and three 
ground-based observations in the J (1.2 gm), H (1.6 um), and Ks (2.1 um) 
bands 14 . Four models fitting the observations are shown in the coloured solid 
curves in main panel, and the coloured circles are the channel-integrated model 
points. The corresponding temperature profiles are shown in the inset. The 
molecular compositions are shown as number ratio with respect to molecular 
hydrogen; all the models have C/O - 1. The thin gray dotted lines show 
blackbody spectra of WASP-1 2b at 2000 K (bottom), 2500 K and 3000 K (top). 
A Kurucz model 30 was used for the stellar spectrum, assuming uniform 
illumination over the planetary disk (i.e weighted by 0.5; Ref 10). The black solid 
lines at the bottom show the photometric bandpasses in arbitrary units. The low 
fluxes at 3.6 and 4.5 gm are explained by methane and CO absorption, 
respectively, required for all fitting models. The high flux in the 5,8 urn channel 
indicates less absorption due to H 2 0. The observations can be explained to 
high precision by models without thermal inversions. Models with strong thermal 
inversions are ruled out by the data (see Figure 3). The red model features a 
thermal inversion at low pressures ( P < 0.01 bar), but the corresponding 
spectrum is almost indistinguishable from the orange model, which has identical 
composition, and identical thermal profile below the 0.01 bar level as the red 



model, but does not have a thermal inversion above the 0.01 bar level. Thus, 
any potential thermal inversion is too weak to be detectable by current 
instruments. 

Figure 2: Constraints on the atmospheric composition of WASP-1 2b. The 
distributions of models fitting the 7 observations (Figure 1) at different levels of 
X 2 are shown. The coloured dots show % 2 surfaces, with each dot representing 
a model realization. The purple, red, green, and blue colours correspond to 
models with / v 2 less than 7, 14, 21, and 28, respectively. Mixing ratios are 
shown as ratios by number with respect to H 2 . At 3a significance, the 
constraints on the composition are H 2 O/H 2 < 8 x 10 “ 6 , CH 4 /H 2 a 4 x 1 0' 6 , CO/H 2 
s 9 x 1 0' 5 , CO 2 /H 2 s 5 x 1 0' 6 , and C/O > 1. The compositions of the best-fitting 
models (with x 2 < 7) span H 2 O/H 2 = 1CT 10 - 1 0' 6 , CO/H 2 = 1 CT 5 - 1 0' 3 , CH 4 /H 2 = 8 
x 1 0‘ 6 - 2 x 10 ' 4 , and C0 2 /H 2 = 5 x 10" 7 - 5 x 10 ' 6 ; the corresponding ranges in 
the elemental abundances are C/O = 1 - 2.5, C/H = 10 3 ~ 1 0 ' 3 and O/H = 10 ‘ 5 - 
IQ ' 3 The solutions with the lowest C/H and O/H ratios correspond directly to the 
lowest CO abundance. Based on thermo-chemical equilibrium, the inferred 
CH 4 /H 2 and H 2 0/H 2 mixing ratios are possible only for C/O a 1, consistent with 
our detection of C/O > 1. The last pane! shows the constraints on the day-night 
energy redistribution, given by (1-A)(1-f r ), where A’ is the bond albedo and 'f r ' is 
the fraction of incident energy redistributed to the night side. Up to f r - 0.45 is 
possible (for A = 0), Thus, the observations support very efficient redistribution. 
An additional observation in the z : (0.9 pm) band was reported recently 15 . 
However, the observation implies a value for the orbital eccentricity inconsistent 
with other data in the literature 13 We therefore decided to exclude this 
observation from the analysis presented here, although including it does not 
affect our conclusions regarding the value of C/O or the temperature structure. 
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Figure 3: Thermal profiles of WASP-1 2b. The purple, red. green, and orange 
profiles correspond to models that fit the observations to within % 2 of 7, 14, 21, 
and 28, respectively (same models as in Figure 2); only 100 randomly chosen 
profiles for each x 2 level are shown, for clarity. The thick, black, solid (dashed) 
curve shows a published profile from a self-consistent model of WASP-1 2b with 
(without) a thermal inversion, adapted from Ref. 20, which assumes solar 
abundances. If a thermal inversion is present in WASP-12b, it is expected to be 
prominent as shown by the solid black curve. A prominent thermal inversion 
between 0.01 - 1 bar is ruled out by the data at 4a. The ostensibly large 
inversions in the figure are at low pressures (below 0.01 bar), which have low 
optical depths, and hence minimal influence on the emergent spectrum (see 
Figure 1). The observations are completely consistent with thermal profiles 
having no inversions. Small thermal inversions are also admissible by the data, 
and could potentially result from dynamics. The thin dotted, dashed, and dash- 
dot lines in black show condensation curves of TiO at solar, 0.1 x solar, and 10 
x solar composition 20 . 
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SI Figure 1: Contribution functions for representative WASP-12b models 
showing the atmospheric origin of flux observed in each bandpass. Two 
representative temperature profiles are shown in panel a (same as the purple 
and green profiles of Figure 1). The contribution functions in panel b (panel c) 
correspond to the solid (dashed) temperature profile in panel a, colour-coded by 
bandpass. The maximum contribution to the emergent flux of WASP- 12b in all 
the channels comes from the lowest layers of the observable atmosphere, 
below the 0.1 bar level. 




